Valley polarization effects on the localization in graphene Landau levels 
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Effects of disorder and valley polarization in graphene are investigated in the quantum Hall 
regime. We find anomalous localization properties for the lowest Landau level (LL), where disorder 
can induce wavefunction derealization (instead of localization), both for white-noise and gaussian- 
correlated disorder. We quantitatively identify the contribution of each sublattice to wavefunction 
amplitudes. Following the valley (sublattice) polarization of states within LLs for increasing disorder 
we show: (i) valley mixing in the lowest LL is the main effect behind the observed anomalous 
localization properties, (ii) the polarization suppression with increasing disorder depends on the 
localization for the white-noise model, while, (iii) the disorder induces a partial polarization in the 
higher Landau levels for both disorder models. 

PACS numbers: 73.43.-f, 73.23.-b, 73.63.-b 



Graphene, a condensed matter stage for quantum elec- 
trodynamics, has attracted an increasing amount of in- 
terest and now monolayer graphite has turned into the 
most versatile electronic system under study PQ. In the 
presence of magnetic field B, an anomalous integer quan- 
tum Hall effect (QHE) [2j [3] is observed in graphene, 
where the sequence of plateaus is shifted by half-integer 
if compared to the usual QHE, as have been theoretically 
predicted [4, 5 . The anomalous sequence of plateaus is 
given by: a xy = (4e 2 /h)(n + 1/2), where integer n is the 
Landau level (LL) index and the factor 4 comes from spin 
and valley degeneracies. 

The unique feature of a Dirac-cone-like band structure 
in the energy range of interest (Fermi energy) has led 
to an avalanche of intriguing observations and a bunch 
of novel predictions 1 j. There are, in particular, two 
important discussions taking place recently in the litera- 
ture on graphene: (i) First, the discussion about the role 
played by valley degree of freedom - started after the ob- 
servation [6] of the lifting of valley degeneracy for the 
lowest (n=0) LL at high magnetic fields. Many theoreti- 
cal works have been studying the mechanisms related to 
valley polarization [7] H M EH EH El E3] and propos- 
ing possible applications [HI E3 EE]- (ii) Second, the 
discussion about anomalous phenomena on the conduc- 
tion and localization around Dirac point. On the ex- 
perimental side, the suppression of weak localization was 
observed in graphene [T7], while on the theoretical side, 
intriguing observations of disorder causing the increase 
(rather than decrease) of the transmission at the Dirac 
point have been reported [TSJ ESI US HI], still without a 
consensus about the origins of this striking mechanism. 

We report here findings in both (i) and (ii) subjects 
of discussion mentioned above and also a connection be- 
tween them, as we found that mixing of valley polarized 
states can explain the anomalous localization we observe 
for graphene systems in the quantum Hall regime. Re- 
lated to this context, although not for graphene systems, 
we call attention to a recent experimental report on the 



influence of valley polarization on the transport proper- 
ties of a two-dimensional electron system [22] . 

The fascinating graphene landscape can be widely ex- 
plored within a single electron framework. The energy 
dispersion relation of graphene is linear around the two 
incquivalent Brillouin zone corners K and K', forming 
two cones in the band structure, whose apex is named 
Dirac-point. We calculate here the density of states 
(DOS), localization properties and the sublattice polar- 
ization of wave functions around Dirac point in the pres- 
ence of a perpendicular magnetic field, within a tight- 
binding model description. The graphene band struc- 
ture is already well described by a nearest neighbor hop- 
ping between s-like orbitals in a honeycomb lattice: H = 

J2i £ i4 C i + t J2 < iJ > ( e ^ ijc l C ]+ e ~ l ' t ' i: ' C j C ^' Wnere C * is tnC 

fermionic operator on site i. The magnetic field is intro- 
duced by means of the phase <pij = 2n(e/h) J\ Adl in the 
hopping parameter t (t^2.7eV for graphene, but we give 
results parameterized by t). In the Landau gauge, tf>ij — 
along the x direction and <fiij= ±7r(a;/a)4>/4>o along the 
TV direction, with <1> /<E»o = Ba 2 V3e/h (a=2.46A is the 
lattice constant for graphene). Here we address two on- 
site disorder models: white-noise fluctuations - emulat- 
ing short-range scatterers - and correlated (smoothed) 
disorder landscapes [23 . For the white-noise model, un- 
correlated orbital energies are sorted within Si < \W/2\, 
while for the correlated disorder, a gaussian correlation 
Si = Ylj £je - ! Ri-Rj !"/ A is assumed, with correlation 
length A and ej < \W/2\. The importance of the disor- 
der in these systems is widely accepted and addressed by 
different groups ES E3 \M HO HD H3 ■ 

The degree of localization of the states is estimated by 
the participation ratio, PR = l/(NJ2iLi IV^I 4 ), where 
ipi is the amplitude of the normalized wavefunction on 
site i, and N is the total number of lattice sites [50] . 
All the results shown here have been calculated on dis- 
ordered unit cells of 60 x 60 sites (hexagonal lattice of 60 
zig-zag chains, each containing 60 sites), with periodic 
boundary conditions. Averages over hundred disorder re- 
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FIG. 1: (color online) (a)(b) DOS for different disorder 
strengths (broadening of LLs increases with disorder) and 
(c)(d) the corresponding participation ratio of states. First 
column (a)(c) corresponds to the white-noise disorder and the 
second (b)(d) to correlated disorder with A=3a. Three LLs 
are shown (n=0,l,2) for flux "l , /$o=l/30. Localization prop- 
erties (observed from the evolution of PR with disorder) are 
very different in the n=0 LL compared to the others. 



alizations are undertaken. For the correlated disorder we 
take A=3a. We consider a low magnetic flux through the 
lattice unit cell, $/<I>q = 1/30, and only few LLs, to keep 
within the continuum limit of the Hofstadter spectrum 
for the hexagonal lattice [3T] . 

Figure 1 shows the DOS and participation ratio (PR) 
around the Dirac point for increasing disorder in two dif- 
ferent models: white-noise (Fig. la and lc), and corre- 
lated disorder (Fig. lb and Id). From the DOS, one can 
see that the energy range shown corresponds, for this 
magnetic flux, to the lowest three LLs (n=0,l,2), well 
resolved at small disorder amplitudes. The spectrum is 
symmetrical around zero energy, so we are not showing 
the hole side (n=-l,-2). The PR, by definition [5U], gives 
the proportion of the N lattice sites over which the wave- 
function is spread. In this way, the smaller the value of 
PR, the more localized the state is. In fact, a character- 
istic localization is observed at LL tails, while extended 
states are identified by peaks in the PR around the center 
of each LL [33]. We first want to call attention in Fig.l 
to the effect of levitation of extended states [29, 32 . In- 
creasing disorder leads to LL broadening (as expected 
due to intra LL mixing) as well as to a Landau band 
repulsion, due to inter LL mixing [331 134] . This band re- 
pulsion can be clearly seen for the DOS of the LL n=l, as 
a shift of the Landau band in direction to E=0. Looking 
at the PRs, one can observe that the delocalized states 
of the 7i=l LL do not move towards lower energies as 
fast as the DOS, characterizing an effective levitation of 
extended states, as for Schroedinger-like carriers at the 



electronic structure extremes [23] . Such levitation is nec- 
essarily frustrated at the n=0 LL. 

The main effect we want to point out in Fig.l is an 
anomalous increase of PR with increasing disorder, oc- 
curring for states of the n=0 LL. This can be observed in 
both disorder models, but is more evident for the white- 
noise case. From Fig.lc, one easily observes PR peaks 
corresponding to the n=l and n—2 LLs having smaller 
values with increasing disorder, as expected, while for the 
n=0 LL, PR around the peak is substantially increased 
for disorder amplitudes from W/t=0.05 to 0.5 (PR values 
at n=0 LL peak are followed as a function of disorder in 
Figures 5 and 6, so more details are discussed further in 
the text). The observed increase in the PR is surprising 
and called anomalous because increasing disorder is ex- 
pected to have the opposite effect: making states more 
localized. The same anomalous effect persists for larger 
system sizes (not shown here). Nevertheless it is worth 
mentioning that size effects would not be able to explain 
this observed increase of PR wi t h increasing disorder. 
In this way, one has to look at another degree of freedom 
in the problem. 

A unit cell of the graphene hexagonal lattice contains 
two carbon atoms, defining two sublattices, A and B. For 
the ideal graphene system (no disorder) in the presence 
of magnetic field, there is a great difference between the 
wavefunctions from the n=0 LL and that from higher 
levels: in the lowest LL, wavefunctions have non-zero 
amplitudes on only one of the sublattices [TO] . But how 
does disorder affect this picture? How does the sublat- 
tice polarization of the n=0 LL evolves as the disorder 
increasingly introduces Landau level coupling? Do the 
higher LLs remain unpolarized with increasing disorder? 
We show next that the role of disorder in destroying sub- 
lattice polarization in the n=0 LL reveals to be related 





FIG. 2: (color online) Wavefunction, and its sublattice de- 
composition, of a valley polarized state from n—0 LL (top) 
and a valley unpolarized state from n=l LL (bottom). Wave- 
function amplitudes (l^l 2 ) are firstly shown over all the 60x60 
lattice sites, and then contributions from each sublattice 
(|^a| 2 and |$s| 2 ) are plotted separately. Results are for a 
correlated disorder with W/t— 1.0, A=3a and <3?/<I>o=l/30. 
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to the anomalous behavior of the PR at n=0 LL. 

We start showing in Figure 2 examples of wavefunc- 
tions probability densities (l^l 2 ) and their decomposi- 
tions in the two sublattices, \*S> a\ 2 and |\I/b| 2 . We de- 
fine here I^E^s)! 2 as the sum of wavef unction amplitudes 
only over the sites of sublattice A(B). In this way, due 
to wavefunction normalization, a\ 2 + l^sl 2 = 1- The 
two decomposed states shown in Fig. 2 correspond to typ- 
ical extended states, one from the n=0 LL (top) and the 
other from the n=l LL (bottom), for a correlated dis- 
order potential. We can clearly see that the state from 
n=0 LL is sublattice (or valley) polarized for the small 
disorder amplitude considered here: the sublattice A con- 
tribution to the chosen state is still very small and can 
be neglected (all plots have the same amplitude scale). 
On the other hand, the state from the n—1 LL is valley 
unpolarized: by eye inspecting the probability densities, 
the contribution of both sublattices are equally relevant. 

The distribution of wavefunction amplitudes over only 
one of the graphene sublattices, or preferentially over one 
of them, is what we call here sublattice (valley) polariza- 
tion. To analyze the limits where this polarization takes 
place and address the questions posed above, we define a 
simple quantity to infer the valley polarization: 



Polarization 



I*/ 



(1) 



This quantity gives the difference in modulus of the con- 
tribution from each sublattice to the total probability 
density of the state. Zero polarization means that the 
state is equally distributed on both sublattices, while to- 
tal polarization, corresponding to values 1.0, means that 
the state is completely on sublattice A or B. 



With this definition, valley polarization of the states 
can be followed for increasing disorder in the white-noise 
and correlated models - the results are shown in Figures 
3 and 4, respectively. Note that the main graphs have 
exactly the same energy range shown in Fig.l. Zoom 
graphs show the surprising fine structure of polarization 
on the n=l and n=2 LLs. Comparing disorder models 
for equivalent LL broadening, it is clear that white-noise 
disorder (Fig. 3) induces more valley mixing on the n=0 
LL than the correlated disorder (Fig. 4). This effect is in 
connection to the fact that the larger the length scale of 
the disorder compared to the lattice constant, the smaller 
the capacity of disorder to cause inter- valley mixing 1 10] . 
However, the wavefunctions polarization observed from 
Figs. 3 and 4 show other intriguing fingerprints that go 
beyond natural expectations, as we discuss next. 

For white-noise model (Fig. 3), the states from the n=0 
band clearly show an interesting " U-shape" polarization 
profile: states at the band tails - strongly localized (as ob- 
served in Fig.lc) - remain more polarized than the delo- 
calizcd ones, from LL center, which become increasingly 
unpolarized with disorder. This relation between local- 
ization length and polarization of the states in the n—0 
LL is not present when the disorder potential is smoother 
(Fig. 4). From Fig. 4 we see that the destruction of the 
sublattice polarization with Gaussian-correlated disorder 
is a slow process, that affects all states from the n=0 LL 
in the same manner, regardless of being localized or de- 
localized. 

Another striking observation that emerges from this 
analysis is the induced polarization beyond the n=0 LL: 
states that are completely unpolarized in the absence of 
disorder, start to show an increasing partial polarization 
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FIG. 3: (color online) Valley polarization of the states within 
n=0,l,2 LLs (same energy range shown in Fig.l), for three dif- 
ferent disorder amplitudes in the white-noise disorder model . 
In the detail, the zoom of the region inside the dashed line, 
corresponding to the fine structure of the polarization in the 
n=l and n—2 LLs. 



FIG. 4: (color online) Valley polarization of the 
states within n=0,l,2 LLs (same energy range shown 
in Fig.l), for three different disorder amplitudes in the 
gaussian-correlated disorder model (A=3a). In the detail, the 
zoom of the region inside the dashed line, corresponding to 
the fine structure of the polarization in the n=l and n=2 LLs. 
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(up to «6%) in the disordered limit. Zoomed graphs 
show that the fine structure of this n > polarization 
presents the " U-shap^ polarization profiles for both dis- 
order models, observed for both n=l and n=2 LLs. 

The results shown in Figs. 3 and 4 are averages over 
hundreds of disorder configurations. However, focusing 
on individual disorder realizations one observes that (not 
shown here) the polarization may jump randomly from 
one sublattice to the other for neighbour states in en- 
ergy within the LL. This behaviour is inferred from the 
quantity described in eq.(l) when the modulus is skiped. 

In Figures 5 and 6 we show in more detail the be- 
havior of the wavefunctions at the center of the n=0 LL 
with increasing disorder, for white-noise and for corre- 
lated disorder, respectively. The graphs show the sub- 
lattice polarization (upper panels) and the participation 
ratio (lower panels) for the delocalized state at E=0 (we 
take in fact mean values of the group of states closest to 
E=0 to minimize the small fluctuations still present even 
when averages over hundreds of disorder realizations are 
undertaken). The anomalous increase of PR with dis- 
order we pointed out in Fig.l is more clearly appreci- 
ated in Figs. 5 and 6. The important point here is that 
the evolution of the PR can be simultaneously followed 
with the valley polarization of these wavefunctions. In 
this way, a correspondence between the degree of local- 
ization and valley (sublattice) polarization can be fully 
appreciated in these figures, as well as the qualitative 
difference between both disorder models. For both mod- 
els, the anomalous increase in the PR is accompanied by 
an equivalent rate of decrease in the polarization. Verti- 
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FIG. 5: (color online) Values of valley polarization and partic- 
ipation ratio (PR) for the states at Dirac point (at the center 
of the n=0 LL) , as a function of white-noise disorder ampli- 
tude. Dashed line indicates the W/t for which the broadening 
of LLs is such that the n=0 LL starts to overlap with neigh- 
bors bands. 



cal dashed lines indicate the disorder amplitude W/t for 
which the n=0 Landau band begin to overlap with the 
n=±l bands (the band tails start to touch each other) 
within each model. For white-noise, Fig. 5, the polariza- 
tion drops faster and is almost suppressed at W/t = 1. 
Just before this threshold in the density of states land- 
scape, the PR reaches its maximum and for higher degree 
of disorder the localization recovers the usual behavior of 
the other Landau levels. On the other hand, for Gaus- 
sian correlated disorder, Fig. 6, the overall picture seems 
to be more involved. The polarization is almost unaf- 
fected at W/t = 2, the Landau band overlap threshold 
for this disorder model. It should be noticed that the 
polarization is still appreciable for disorder strengths at 
which any modulation in the density of states are com- 
pletely washed out. Thi s is the case for W/t > 3.5 (see 
Fig. lb), but the polarization at W/t = 4.0 is not negligi- 
ble (around 65 %). The observation that the PR contin- 
ues to increase up to these high disorder values seems to 
be related to the fact that the valley polarization is still 
appreciable for this disorder. 

In conclusion, these results represent a strong evidence 
of connection between the anomalous localization proper- 
ties observed for the n=0 LL and the way that disorder 
frustrates the valley polarization in this Landau level. 
As disorder increases in the graphene system, sublattice- 
polarized wavefunctions from the n=0 LL increasingly 
spread to the other sublattice, allowing more sites to 
participate in the wavefunction. This increases, by defi- 
nition, the participation ratio of the state. The expected 



Gaussian-Correlated Disorder 




FIG. 6: (color online) Values of valley polarization and partic- 
ipation ratio (PR) for the states at Dirac point (at the center 
of the n=0 LL), as a function of Gaussian-correlated disor- 
der amplitude. Dashed line indicates the W/t for which the 
broadening of LLs is such that the n=0 LL starts to overlap 
with neighbors bands. 
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and well known localization with increasing disorder only 
comes back to this scenario when disorder has signif- 
icantly destroyed the sublattice polarization. Another 
striking found is that the polarization suppression, as a 
function of energy within the lowest LL, is model depen- 
dent, with the intriguing "U-shape" for the white-noise 
model. On the other hand, the increasing LL mixing also 
induces a partial polarization in the higher levels. These 
observations may establish a key observable for unveil- 
ing the striking differences between the properties of the 
Dirac-like in comparison to Schroedinger-like QHE. 

We are grateful to A.H. Castro Neto for stimulating 
discussions and L. Brey for calling the attention to the 
problem. This work is supported by FAPESP and CNPq. 
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